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Abstract 



Two points on the Coulomb branch of M = 4 super Yang Mills are investigated using their supergravity 
duals. By switching on condensates for the scalars in the N = 4 multiplet with a form which preserves 
a subgroup of the original i?-symmetry, disk and sphere configurations of D3-branes are formed in the 
dual supergravity background. The analytic, canonical metric for these geometries is formulated and 
the singularity structure is studied. Quarks are introduced into the corresponding field theories using 
D7-brane probes and the meson spectrum is calculated. For one of the condensate configurations, a 
mass gap is found and shown analytically to be present in the massless limit. It is also found that 
there is a stepped spectrum with eigenstate degeneracy in the limit of small quark masses and this 
result is shown analytically. In the second, similar deformation it is necessary to understand the full 
D3-D7 brane interaction to study the limit of small quark masses. For quark masses larger than the 
condensate scale the spectrum is calculated and shown to be discrete as expected. 
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1 Introduction 



Substantial effort has gone into understanding the properties of gauge theories dual to supergravity 
backgrounds which asymptote to AdSs x S 5 . This progress has been made possible by the conjectured 
AdS/CFT correspondence [H[H[3]. In particular, there has been much interest in attempting to obtain 
QCD-like models by breaking the supersymmetry and conformal symmetry via the inclusion of relevant 
operators [U O [6j [7] . The addition of quarks [H ] has also heralded a great leap in creating realistic 
models and has allowed us to calculate many non-perturbative quantities. Over the last year several 
toy models of five-dimensional holography have also made progress in describing theories with a small 
number of colours and even the most naive of these scenarios appears to give remarkable agreement with 
lattice QCD and experimental observations of meson masses and decay constants [271 128, 29, 30], I3H [32] . 
Generally, even in the simplest deformations, calculations must be performed numerically in order to 
find solutions to five and ten-dimensional equations of motion. Similar models of AdS slices were also 
considered in [331 US] where glueball spectra were studied. 

In this paper we study one particular deformation [35] of the AdSs x S 5 geometry using D7-brane 
probes. The field theory dual to this geometry retains the full N = 4 supersymmetry but breaks 
the SU{4)r symmetry by the addition of condensates for the three complex scalar fields. Two point 
correlation functions and Wilson loops have been studied in this theory on the Coloumb branch of 
N = 4 Super Yang-Mills and the features of the scalar spectrum and screening have been explained 
in terms of ensembles of brane distributions. In the current work, the supersymmetry allows the 
analytic form of the metric to be obtained which encodes the field theory in its canonically normalised 
form, a task which is more difficult in non-super symmetric backgrounds. In the canonically normalised 
coordinates the meson spectrum can be calculated as a function of quark mass once D7-branes have 
been introduced. 

The addition of fundamental matter to the theory under consideration has been studied briefly 
in [24] . When the D7-brane flow was calculated in this geometry, there appeared to be a small 
but finite quark bilinear condensate for non-zero quark mass. This result is surprising because this 
supersymmetric theory should not support a chiral condensate. This anomalous result is however due 
to the use of an unsuitable basis with which to describe the geometry in a holographic setting. The 
details of this are addressed in section 14.11 It is trivial to prove that any supersymmetric background 
when written in canonical form will define a stable field theory with zero vacuum expectation value 
for the quark bilinears. It was shown using this background as a simple example that the original, 
geometric interpretation of chiral symmetry breaking was not sufficient for the analysis of backgrounds 
out of their canonical form. A method was developed by which the potential felt by a D7-brane in the 
singular region could be studied out of canonical form. It is the aim of the current paper to find the 
canonical basis with which to describe the field theory and study the spectrum therein. 

2 The Supergravity Geometries 

The two geometries of interest in this paper are two of a set of five solutions discussed in [35] which are 
all asymptotically AdSs x S 5 and are sourced purely by D3-branes. Each background is formulated in 
terms of a D3-brane density distribution function. 

It is possible to find the analytic, canonical form for all five of these supergravity backgrounds 
though only in two of the cases is the form of the metric simple enough to calculate the meson 
spectrum. 

Each background is dual to an M = 4 field theory with a scalar condensate, preserving a subgroup 
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of the original SU(4)r symmetry. In each case the metric is given by 



ds< 



H 



dx 



M 



(1) 



where the warp factor is 



H 



t 

Jo 



L 4 

d n ua(u})— — 

l<3|<i \y-vr 



(2) 



u is a vector in n of the six dimensions transverse to the D3-brane worldvolume and y is a vector 
in all six of these dimensions. I parametrises the size of the D3-brane distribution and is the single, 
extra, free parameter in each of these geometries. The integral is performed over the region of space 
with support from the distribution function. The dimension of the distribution, n, together with the 
density function, the preserved symmetry and the form of the scalar condensate are provided in table 

m 



n 




Preserved Symmetry 


Scalar Condensate 


1 
2 
3 
4 

5 


l l 

n 2 l 2 

i ( 5{i 2 -uj 2 ) e(i 2 -u; 2 ) \ 


50(5) 
50(4) x SO{2) 
50(3) x 50(3) 
50(2) x 50(4) 

50(5) 


^(1,1,1,1,1,-5) 
i(l,l,l,l,-2,-2) 
±(1,1,1,-1,-1,-1) 
1(2,2,-1,-1,-1-1) 

^(5,-1,-1,-1,-1,-1) 



Table 1: D3-brane density distribution functions preserving subgroups of the 50(6) symmetry. 
The condensate is a vector in the six-dimensional space of real scalars of M = 4 SYM. 



3 Obtaining the Canonical Coordinates 

The two backgrounds of interest are those which preserve SO (2) x 50(4) subgroups of the original 
50(6) symmetry. These are particularly interesting because once fundamental matter is added, the 
chiral symmetry is described explicitely by the 50(2) geometrical symmetry. 

A non-canonical analytic form for these backgrounds has been given in From this form of 

the metric it is simple to hnd the canonical system. A D3-brane is introduced which describes a field 
theory with six scalar fields. The canonical metric is defined as that in which the six scalar fields are 
simultaneously canonically normalised 0. 

Note that recent work [39] describes an unambiguous method of finding the natural coordinate 
system for supersymmetric deformations using holographic renormalisation. 

One of the backgrounds of interest dual to M = 4 SYM with an adjoint scalar condensate has 
been studied previously \A2\ [2"3] in the context of flavoured holographic models. This geometry was 
conjectured to be equivalent to the n = 2 geometry of tabled! though we show in this section that it 
is really another parametrisation of the n = 4 deformation. 

1 Thanks to K Sfctsos for pointing out that this canonical coordinate system has been calculated previously 
in [37] and [35] 
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The original coordinate system used to study mesons in this background is a limiting case of the 
J\f = 2 background of [6]. To return to the full J\f = 4 theory, the two five-dimensional supergravity 
scalar fields are equated to acquire a theory with six scalar vevs of the form in table HJ In this particular, 
limiting case of the supergravity solution, the dilaton becomes constant. It may be interesting to study 
those geometries with a running dilaton in their canonical coordinates in the future. 

In this supergravity background, there are two fields of interest. One is the five-dimensional scalar 
field, X-, an d the other is the warp factor, A, multiplying the Minkowski space-time components in 
the five-dimensional truncation of the ten-dimensional metric. In [35] the lift of the five-dimensional 
supergravity theory was obtained and we use the resulting metric in what follows. 

The equations of motion for the scalar field and warp factor are 



and 

2A z a, 

B? y 6 - 1 ■ { 1 



I 2 X 4 



R 2 x 6 

The metric for this background in the unphysical coordinates is given by 



ds l = ^-e ZA dx z M - ^— \du z + U6 2 + -j^dp + dftij 1 , (5) 

where 

X = cos 2 6 + x 6 sin 2 6 . (6) 

By probing with a D3-brane, the action for the six scalar fields is seen not to be canonically normalised, 
though the moduli space is manifest. In this case it is possible to use the first order supergravity 
equations of motion to find the correct coordinate system in which to describe the field theory in its 
canonical form. The result of this transformation is given by 

ds 2 = H~Ux 2 M - m {dp 2 + dull + dul + p 2 dnl) , (7) 

where the warp factor is 

H(p, u, I) = — . (8) 

Z 4 + 21 2 (lo 2 - p 2 ) + (a; 2 + p 2 ) 2 + {I 2 + lo 2 + p 2 ) ^(to 2 + (I - p) 2 ){u 2 + (I + p) 2 ) 

H is plotted in the (u, p) plane in figure[T]in order to illuminate the singularity structure. This deforma- 
tion corresponds to an S* 3 distribution of D3-branes spanning the locus p = i in the R 4 = S 3 x p plane. 
This is therefore the analytic, canonical form of the n = 4 metric in [35] which describes an M = 4 field 
theory with vacuum expectation values for all six scalar fields in the configuration (2, 2, —1, —1, —1, —1). 
The original metric also encodes the n = 2 solution with scalar vev (1, 1, 1, 1, —2, —2), however this 
configuration has negative I 2 as discussed in [35] . 

The warp factor for the n = 2 configuration can be plotted in the same way and is illustrated 
in figure [2j This corresponds to a D 2 distribution spanning the (ws,^) two-plane. It is possible 
to calculate the canonical, analytic form for the other three metrics in [35] which preserve different 
subgroups of the original SO (6) symmetry. However, the analytic form for these backgrounds are too 
complicated to calculate the meson spectra so will be discussed no further. 
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Figure 1: Singularity structure of the warp factor for the n = 4 geometry in the (u),p) plane 
with I 2 = -1. 




Figure 2: Singularity structure of the warp factor for the n = 2 geometry in the (u, p) plane as 
in figured] but with I 2 — 1. 



4 Mesons from D7-Brane Probes 

Having found the correct coordinate system in which to describe a canonically normalised field theory 
living on the D3-branes, we can study the theory with the addition of quarks. 
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We start with the background described by equation [7] with positive I 2 , which preserves an SO (4) x 
50(2) subgroup of the original symmetry. The D7-brane is embedded by filling the four Minkowski 
space-time directions, the p-direction and wrapping the S 3 . This wrapped cycle ensures the stability 
of the brane configuration. 

Before calculating the ^-dependent excitations of the brane, we must study its flow purely in the 
p-direction. This corresponds to calculating oj§ and ujq as a function of p. However, because of the 
super symmetric nature of this background, the warp factor in front of the R 6 is the inverse of that in 
front of the Ai 4 plane. This means that the x-independent flow is exactly the same as in the AdSs x S 5 
background to which the solution is known analytically. 

We can solve the equation of motion and obtain the following solutions. 

J vp -4 



where u> = \J to 2 + uj\ . It is clear that this function will not be real all the way to p = for c\ ^ 0. 
The physical solution, corresponding to the renormalisation group flow of the brane, must therefore 
be oj = C2, equivalent to a quark mass but no quark bilinear condensate. This will always be the case 
for a supersymmetric solution where the warp factors cancel leaving the AdSs x S 5 equation of motion 
(excluding ^-dependent fluctuations). 

Because of the manifest 50(2) symmetry between lo§ and loq, we are free to choose the direction in 
which to explicitly break this rotational invariance. For simplicity, we choose the solution uj^ = m and 
u>6 = 0. We now want to study the mesonic fluctuations about this brane flow. We study the modes in 
the a;6-direction given by ujq = + Cjq(p, x). Note that because there is no chiral symmetry breaking 
in this background, the mesonic excitations in the L05 and o^-directions will be identical. This means 
that the positive and negative parity states will be degenerate. 

The action for ujq up to quadratic order is given by 



S 



<9^6 X 2 
dp 



+ 



d£ 6 \ 2 2R 4 



Ox J /4 + 2 /2( m 2 _ p 2) + ( m 2 + p 2)2 + (J2 + m 2 + p 2)^( m 2 + (/ _ p)2)( m 2 + (/ + p )2) 



(10) 



For small oscillations about the flow us(p), the meson interaction terms will be subdominant and the 
function ojq can be treated as a plane wave in the Minkowski space-time directions and therefore the 
ansatz for this function is given by 

^6 = Me*" . (11) 

This ansatz which is independent of the three-sphere coordinates corresponds to an i?-singlet, spinless 
meson wavefunction. The equation of motion for f{p) is given by 

2M2 ^ 4/( "' , + 3/-M + p/"M = , 

Z 4 + 2l 2 (m 2 - p 2 ) + (m 2 + p 2 ) 2 + {I 2 +m 2 + p 2 )^/(m 2 + (I - p) 2 )(m 2 + (I + p) 2 ) 

(12) 

where M 2 = —k 2 . The eigenvalues, M, are given by the values for which the flow of / is well behaved 
all the way to p = and normalisable in the UV. It must also have the correct scaling dimensions in the 
AdSs x S 5 limit to describe a mesonic excitation. Of the two UV solutions, the solution corresponding 
to meson fluctuations, as opposed to an x-dependent mass is f(p) — ► -%. 
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It is possible to remove the explicit R and I dependence in the above equation of motion by 
performing the following rescaling: 

m = fhl, p = pi, MR 2 = Ml, (13) 
giving us the following equation of motion 

1 + 2(m 2 - p 2 ) + (m 2 + p 2 ) 2 + (1 + m 2 + p 2 )^(rh 2 + (1 - p) 2 ){m 2 + (1 + p) 2 ) 

(14) 

Before performing this rescaling we can take the I — > limit and find that the numerical values 
coincide with the known analytic results of the pure AdSs x S 5 spectrum }41j : 

2m 



M = -jpV ( n + l)(n + 2) , n = 0,l,... (15) 

showing that in this limit the numerics are under control. For the rescaled equation (eq !14p we can 
study the lowest mass meson as a function of the quark mass. Figure [3] shows the value of the first 
meson mass as a function of the quark mass. The important point to note here is that there appears 
to be a mass gap in the m — ► limit for I ^ 0. The numerics make this calculation difficult, though 
at m = 10~ 10 the value of Mi is 0.28. Note that in contrast to the equation of motion with I = 0, 
the D7-brane equation is perfectly well behaved in this limit and has discrete eigenvalues. This will be 
shown analytically in section [4721 




Figure 3: Mi against log 10 m for I = 1,R = 1 (solid line). This becomes exactly degenerate 
with Mi for I = 0, R = 1 (dashed line) in the large m limit. Note the appearance of a mass gap. 
The scale is set by the AdS radius R which can be tuned by hand to compare with lattice data. 



Having calculated the first meson mass, we can study the spectrum for the higher excited states. 
The results of this are shown in figure [H We can see from figure H] that there is a degeneracy for 
quark masses much less than the condensate scale (given by the gap between the singularities in the 
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0.5 



m 



1 



1.5 



Figure 4: M n against m for the first seven states. For large quark masses the behaviour is 
simple as the brane cannot resolve the double singularity structure, in the small mass limit 
there is a degeneracy of states. In the large quark mass limit the spectrum goes like MR 2 = 
2m-\/ (n + l)(n + 2) just as in the pure AdS 5 x S 5 case. 



supergravity geometry). For large quark masses, corresponding to the D7-brane lying far from the 
singular region, there is a spectrum which is given by 

MR 2 = 2m^(n + l)(n + 2), n = 0, 1, ... (16) 

as expected because the brane cannot resolve the separated singularities for large quark masses. For 
smaller masses however we see that pairs of excited states become degenerate and indeed there appears 
to be an exact degeneracy as we get to the massless limit. The lowest state however does not appear 
to have a degeneracy and we will show in section 14.21 that this mode does not appear in the exact 
massless limit when we study the equations analytically. 

4.1 Significance of the Canonical Coordinates 

It is important at this stage to note why we are interested in finding one particular coordinate system 
for this problem. Though it appears strange that coordinates matter in what is essentially an eigen- 
value problem we see in the following that this canonical coordinate system simplifies the calculations 
considerably. Without the change of coordinates not only would we be unable to obtain an analytic 
expression for the meson spectrum in the massless quark limit but because of the numerical behaviour 
of the supergravity field, Xi we would not be able to obtain a reliable spectrum at all. This means that 
the numerical instabilities are critical to the eigenvalue problem. 

These instabilities stem from two regions in the geometry. The first of these is the singular region 
where the supergravity field, asymptotes to infinity. The second region is where the space returns 
to pure AdSs x S 5 . In this region the supergravity field asymptotes to unity but the further we go into 
the ultraviolet, the more significant becomes the accuracy of the solution in insuring a numerically 
stable geometry. This necessity for ever higher numerical accuracy in the ultraviolet can be seen in 
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equations and [H In the canonical basis the singular behaviour cancels explicitly in the flow equation 
for the D7-brane. 

The first hint that these instabilities are critical is seen when we calculate the stable value of 
the quark bilinear condensate as a function of the quark mass. The stable flows of the D7-brane for 
various quark masses and the quark bilinear condensate versus quark mass are plotted in figures [5] and 
[6] respectively. In the canonical coordinates we see analytically that there is no condensate present 




Figure 5: Stable D7-brane flows for varying quark masses. For ease of visualisation the flows 
are plotted using the coordinates r 2 + v 2 = R 2 e 2u and - = tan 8 with a singularity at radius R. 




Figure 6: Quark bilinear condensate versus quark mass calculated from the stable D7-brane 
flows in figure El 



for any quark mass. However, in the original coordinate system (equation [5]) , we see that for finite 
r the stable D7-brane solution does not lie flat with respect to the r axis indicating a condensate 
present for non-zero quark mass. In theory we should be calculating the condensate at infinite energy, 
where the space returns to pure AdSs x S 5 . This is not possible as we have had to solve the equation 
numerically. Solving for finite r will always give a non-zero condensate for non-zero quark mass. As 
we perform the calculation further into the UV where this calculation should be more accurate, the 
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numerical instabilities of the scalar supergravity field become critical and a meaningful result becomes 
harder to calculate. 

When calculating the meson spectrum as a function of the quark mass, the spectrum which is 
obtained is also found to be a chaotic one, clearly influenced by the singular behaviour in the IR 
and the numerical instabilities in the UV, especially in the region of small quark mass which is the 
domain of most interested. This instability is critical to the calculation and means that the interesting 
properties which can be found analytically and are discussed in the next section are not observable in 
the original basis using the computational techniques employed here. 

4.2 Analytical Results 

Though we know from the m — > limit of equation [15] (The AdSs x S 5 meson spectrum) that at m = 
the spectrum becomes continuous, this appears not to be the case from the numerical calculation in 
the I deformation. In this case, we can take the m = limit explicitly in the equation of motion 
and retain a discrete spectrum. 

M 2 pR A f{p) + l\l 2 - p 2 )(3/'(p) + pf"(p)) = . (17) 

The solution to this equation is given in [351 EE] for the case of the glueball spectrum and is given by: 

MR 2 = 2l v / (n + l)(n + 2) n = 0, 1, ... (18) 

(see also [38] for higher angular momentum eigenvalues). It is also precisely the same spectrum as the 
AdSs x S 5 spectrum with an interchange of the quark mass, m and the deformation parameter, I (the 
magnitude of the scalar condensate). It is particularly interesting that at this point on the moduli 
space the spectrum of M = 2 quarks has the same spectrum of states as the M = 4 theory. 

There is an apparent discrepancy between the analytical and numerical values calculated in the m — > 
limit. Interestingly, the masses using the numerical and analytic methods are exactly equal except 
for the very first state. However, because the infinitesimal and exact massless limits are qualitatively 
different we may expect a discrete change in behaviour. 

For the numerical calculation, in the case, m = 1CP 10 and I = R = 1, the spectrum is given by 
the values in table [2[ The degeneracy is given by the number of mesons with approximately the same 
eigenvalue. 



M 


Degeneracy 


0.28 


1 


2.8 


2 


4.9 


2 


6.9 


2 


8.9 


2 



Table 2: The mass spectrum for I = R = 1 and m — > calculated numerically. The degeneracy 
is given by the number of states with converging eigenvalues. 

This should be compared with the exact m = result obtained analytically where the degeneracy 
is not obtainable. These results are given in table 
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JVI 


Degeneracy 


9 8 


n.a 


4.9 


n.a 


6.9 


n.a 


8.9 


n.a 



Table 3: The mass spectrum for I = R = 1 and m = calculated analytically. The degeneracy 
is not possible to study in the analytic massless limit because it is not possible to calculate the 
number of eigenstates of a given eigenvalue. 



4.3 The n = 2 Spectrum 

As discussed in section [3] this deformation is obtained from the n = 4 deformation with the interchange 
I 2 — > —I 2 . This transformation however has a significant effect on the nature of this geometry from the 
point of view of a D7-brane probe. Again, we can see from the simple product form of the geometry 
with inverse warp factors for the A4 and R 6 that to zeroth order in mesonic excitations, the D7- 
branes will he flat and not notice the singular structure. In figure [21 it was shown that there is a disk 
singularity lying in the (u^,Uq) plane of radius I. This means that a brane corresponding to adding 
quarks of mass m for m < I will pass straight through the singular region. We know that this is not 
a physical configuration and the full interaction between the D3-brane stack and the D7-brane probe 
would be needed to understand this case fully. It is however possible to study quarks with m > i in the 
supergravity, probe approximation. We find that as expected, for large masses, the spectrum returns 
to the AdSs x S 5 values, just as in the large mass limit of the n = 4 geometry. The spectrum down 
to m = I is provided in figure [7] Though it appears from this diagram that the states may become 




1.02 1.04 1.06 



1.08 



1.1 



in 



Figure 7: M n against m for the first seven states of the n = 2 deformation for m > I (solid 
lines) . 



degenerate (and possibly massless) at m = I this is not apparently the case from these numerical 



II 



computations, The first seven states are given by the following values at m = I: 



M n = {1.46, 1.59, 1.77, 2.00, 2.26, 2.54, 2.84} (19) 

In fact the equations are perfectly well behaved for m < I and the spectrum is calculable in this region, 
however the results should not be trusted as clearly interactions between the D7-brane and the singular 
D3-brane distribution must be understood fully away from the supergravity limit. 

We find that M n is linear in I for m = I, though the n-dependence is not as simple as the massless 
n = 4 geometry. Although the equation of motion simplifies significantly to 

= + 3pf'[p]+p 2 f"[p], (20) 



4 + p + p^/l + p 

no analytic solution to this equation was obtained. 

The qualitative difference between the meson spectrum in this geometry and the n = 4 case is that 
there does not appear to be a degenerate spectrum here. The qualitative difference from the point of 
view of the D7-brane is that there is a continuous distribution of branes rather than a set of singular 
points with a separation of 21. 



5 Conclusions 

We have found an analytic, canonical form for two scalar deformations formulated in [35]. The sin- 
gularity structures are exactly as expected from the D3-brane density distribution functions. In these 
two cases the 50(6) symmetry is broken to SO (4) x SO (2) which encodes the chiral symmetry of 
the field theory explicitly, it is possible to calculate the meson spectrum from excitations of probe 
D7-branes. Interestingly, in one of these backgrounds, it is possible to find the spectrum analytically 
in the limit of exactly massless quarks. This result includes the existence of a mass gap proportional 
to the deformation parameter with exactly the same spectrum as the adjoint scalar two point Greens 
function of the M = 4 theory of [35]. A degeneracy is discovered in the limit of small quark masses 
(c.f the deformation parameter I). 

In the second SO(4) x SO (2) preserving background, it is possible to find the spectra for quarks with 
larger masses, m, than the deformation parameter, I. The spectrum in this case does not appear to be 
degenerate though to fully understand the m ~ I limit, higher order corrections to the supergravity limit 
must be calculated. The overall conclusion is that by finding the correct analytic, physical coordinates 
in which to describe the dual gravity theory, we can study the elaborate structure of meson spectra 
in field theories with complicated condensate terms switched on. If we can do the same thing in the 
non-supersymmetric analogues of these geometries, we may be able to gain some more insight into real 
world hadron spectra. 
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